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Abstract 

Elastic and proton-dissociative photoproduction (7^ p^p, p^N, respec- 

tively, with p^ — > 7r+7r~) has been studied in ep interactions at HERA for photon-proton 
centre-of-mass energies in the range 50 < < 100 GeV and for \t\ < 0.5 GeV^, where t is 
the square of the four-momentum transfer at the proton vertex; the results on the proton- 
dissociative reaction are presented for masses of the dissociated proton system in the range 
< O.IVF^. For the elastic process, the Tr+vr" invariant mass spectrum has been inves- 
tigated as a function of t. As in fixed target experiments, the p^ resonance shape is asym- 
metric; this asymmetry decreases with increasing as expected in models in which the 
asymmetry is ascribed to the interference of resonant and non-resonant vr+vr^ production. 
The cross section has been studied as a function of W; a fit to the resonant part with the 
form gives a = 0.16 ±0.06 (stat.) Iq 15 (syst.). The resonant part of the jp — > TT~^Tr~p 
cross section is 11.2 it 0.1 (stat.) '^I'l (syst.) ph at (W) = 71.7 GeV. The t dependence 
of the cross section can be described by a function of the type Apexp (— + Cpt^) with 
bp = 10.9 ± 0.3 (stat.) IJ;^ (syst.) GeV-^ and Cp = 2.7 ± 0.9 (stat.) Ij;? (syst.) GeV"! 
The t dependence has also been studied as a function of W and a value of the slope 
of the pomeron trajectory a^' = 0.23 ± 0.15 (stat.) Igio? (syst.) GeV~^ has been 
deduced. The p*^ spin density matrix elements Tqq, rf^i and 3?e[r5'o] have been mea- 
sured and found to be consistent with expectations based on s-channel helicity conser- 
vation. For proton-dissociative tt~^tt~ photoproduction in the p^ mass range, the distri- 
butions of the two-pion invariant mass, W and the polar and azimuthal angles of the 
pions in the helicity frame are the same within errors as those for the elastic process. 
The t distribution has been fitted to an exponential function with a slope parameter 
5.8 ± 0.3 (stat.) lb 0.5 (syst.) GeV^^. The ratio of the elastic to proton-dissociative p^ 
photoproduction cross section is 2.0 ± 0.2 (stat.) it 0.7 (syst.). 
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1 Introduction 



Elastic photoproduction of p° mesons, 7p — > p^p, has been studied in fixed target experiments 
at photon-proton centre-of-mass energies W up to 20 GeV [l|]-p5l and at HERA ||2^-||2^ for W 



up to approximately 200 GeV. In both cases the reaction exhibits the features of a soft diffrac- 
tive process, namely a weak energy dependence and a differential cross section da/dt oc exp {ht) 
at low \t\ values, where t is the squared four-momentum exchanged between the photon and the 
proton. These features, typical of elastic hadron-hadron interactions, are consistent with the 
expectations of the Vector Meson Dominance model (VDM) pOf in which the photon is assumed 



to fluctuate into a vector meson before scattering from the proton. In elastic p° photoproduc- 
tion the photon thus appears to behave like an ordinary hadron which interacts elastically with 
the proton. Many aspects of p° photoproduction remain however to be clarified - among them, 
the W dependence of the cross section, the origin of the asymmetric p° resonance shape and 
the extent to which the helicity of the photon is transferred to the vector meson. Perturbative 
QCD calculations which have been able to succesfuUy describe the photoproduction of J/ip 
mesons are not strictly applicable to photoproduction at low \t\. In general, photopro- 
duction of p° mesons at HERA may offer a means of investigating the nature of soft hadronic 
interactions as well as the hadronic features of the photon. 

Little is known about proton-dissociative p° production with real photons, 7p P^N, 
where is a state of mass into which the proton diffractively dissociates. Data exist 
for the virtual photon case: the HI Collaboration at HERA has recently investigated proton- 
dissociative p° production for photon virtualities > 7 GeV^ . The HI results indicate that 
the cross section for this process has the same dependence on and W and the same helicity 
structure as the elastic reaction. These observations support the hypothesis of factorisation of 
the diffractive vertex which has been extensively studied in hadron-hadron reactions (see 



e.g. |^-||3^). The HI data also show that the t distribution is exponential but shallower than 
that for the elastic case. In photoproduction, proton-dissociative production of p° mesons can 
provide yet another way to study diffraction, the hadronic properties of the photon and the 
nature of soft hadronic processes. In conjunction with the elastic reaction, it can provide a 
test of factorisation. Moreover, a detailed understanding of proton-diffractive dissociation is 
mandatory for the study of the elastic reaction, for which it is the main source of background 
when the scattered proton is not measured. 

This paper describes a measurement of p" photoproduction in the elastic and proton-dissocia- 
tive reactions. The measurement was performed using data collected in 1994 by the ZEUS 
experiment at HERA for the processes ep evr+vr^p and ep — > e7r~^7f~N at small photon 
virtualities, < 4 GeV^. The symbol e indicates positrons. 

For the data presented here the scattered positron was not detected. The scattered proton 
was measured only for a subsample of the data. In general the relevant kinematic quantities 
were determined from the measured three-momenta of the two pions from the p° decay. 



With respect to the ZEUS 1993 data p6[, the present results on elastic p° photoproduction 



feature larger statistics, a wider W range and smaller systematic uncertainties. Two of the 
main contributions to the uncertainties of the 1993 results were significantly reduced: the 
calorimeter trigger efficiency was evaluated directly from the data and the contamination of 



1 



proton-dissociative events was determined by using a subsample of the data in which elastic 
events were unambiguously selected by detecting the final state proton. 

The larger statistics and wider kinematic range allowed the study of the vr+vr" mass spectrum 
as a function of t, W and the decay pions' polar and azimuthal angles in the helicity frame. 
The shape of the vr+vr" mass spectrum in the reaction ep — > e7T~^7r~p is interesting since, in the 
framework of the Soding model [Q, it depends on the interference between resonant and 



non-resonant tt+tt production. 

The cross section for elastic p° photoproduction, a-yp^pOp, was extracted as a function of W. 



The W dependence of the cross section, in Regge theory ^9[, is related to the intercept ajp(O) 
of the pomeron trajectory exchanged between the proton and the hadronic fiuctuation of the 
photon. 

The differential cross section da/d\t\ was determined and its shape studied as a function of 
W. Regge theory predicts that the slope of the exponential t distribution becomes increasingly 
steep with increasing W; the rate of change of the t slope with W, at high values of W, is 
related to the slope of the pomeron trajectory, a^'. 

The decay pion angular distributions in the helicity frame were studied and the p° spin 
density matrix elements Tqq, r^"^^ and J^efr^Q] determined. The behaviour of these matrix 
elements as a function of the two-pion invariant mass M^r^ and of W was investigated. This 
made it possible to test the validity of s-channel helicity conservation (SCHC). 

We also present results on p° photoproduction with diffractive dissociation of the proton 
in the range < O.IVT^. The limit was chosen following refs. |]32| , ^ and corresponds to 



the region where diffractive interactions dominate. The distributions of M^^^, W, the polar and 
azimuthal angles of the pions in the helicity frame and the t dependence of proton-dissociative p° 
photoproduction were studied. The ratio of the cross sections for elastic and proton-dissociative 
photoproduction was determined. 

Finally, we used the data on the reaction ep eTi^Ti~p to evaluate the pion-proton total cross 
section in a model dependent way. As mentioned earlier, the shape of the tt+tt" mass spectrum 
is sensitive to the interference between resonant p° and non-resonant 7r"'"7r~ production. In the 
latter case, one (or both) of the pions interacts with the proton. The mass spectrum thus 
depends on the pion-proton cross section. In the framework of a calculation by Ryskin and 
Shabelski |^T[], we determined this cross section at a pion-proton centre-of-mass energy of the 



order of 50 GeV, beyond the reach of any existing pion beam. 



2 Experimental set-up 
2.1 HERA 

The data discussed here were collected in 1994 using the HERA collider which operated with 
820 GeV protons and 27.5 GeV positrons. The proton and positron beams each contained 
153 colliding bunches, together with 17 additional unpaired proton and 15 unpaired positron 
bunches. These additional bunches were used for background studies. 
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2.2 The ZEUS detector 



A detailed description of the ZEUS detector can be found elsewhere The components 

which are most relevant for this analysis are briefly discussed below. 

Charged particles are tracked by the inner tracking detectors which operate in a magnetic 
field of 1.43 T provided by a thin superconducting solenoid. Immediately surrounding the 
beam pipe is the vertex detector (VXD), a drift chamber with 120 radial cells, each with 12 
sense wires [^. It is surrounded by the central tracking detector (CTD), which consists of 72 
cylindrical drift chamber layers, organised into 9 superlayers covering the polar angle region^ 



15° < 9 < 164° |H5[. The transverse momentum resolution for tracks traversing all superlayers 



is a{pT)/pT ^ ^(O.OOSpt)^ + (0.016)2, with pt in GeV. The Rear Tracking Detector (RTD) 
consists of a planar drift chamber with three layers of drift cells with the wires oriented at 0°, 
+60° and —60° with respect to the horizontal plane; polar angles between 160° and 170° are 



covered 46 



The high resolution uranium-scintillator calorimeter (CAL) |^ consists of three parts: the 
forward (FCAL), the barrel (BCAL) and the rear calorimeter (RCAL); they cover the polar 
angle regions 2.6° to 36.7°, 36.7° to 129.1°, and 129.1° to 176.2°, respectively. Each part is 
subdivided transversely into towers. The towers are segmented longitudinally into one electro- 
magnetic section (EMC) and one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAC). These sections are further subdivided into cells; each cell is viewed by two photo- 
multiplier tubes. The CAL energy resolution, as measured under test beam conditions, is 
aE/E = 0.18//E for electrons and as/E = 0.35//E for hadrons {E in GeV). 

The Veto Wall, the C5 counter and the small angle rear tracking detector (SRTD) 
all consist of scintillation counters and are located a.t Z = —730 cm, Z = —315 cm and 
Z = —150 cm, respectively. Particles which are generated by interactions of protons with 
residual gas molecules in the beam pipe (proton "beam-gas" events) upstream of the nominal 
ep interaction point reach the RCAL, the Veto Wall, the SRTD and C5 at different times than 
particles originating from the nominal ep interaction point. Proton beam-gas events are thus 
rejected by timing measurements in these detectors. 



The proton remnant tagger (PRTl) |^ is used to tag events in which the proton diffractively 
dissociates. It consists of two layers of scintillation counters perpendicular to the beam and 
is positioned at Z = 515 cm. The two layers are separated by a 1 mm thick lead absorber. 
Each layer is split vertically into two halves and each half is read out by a photomultiplier 
tube. The counters have an active area of 30 cm x 26 cm with a hole of 6.0 cm x 4.5 cm at the 
centre to accommodate the HERA beam pipe. The PRTl covers the range in pseudorapidity 
(7] = - In tan (9/2)) from 4.3 to 5.8. 



The Leading Proton Spectrometer (LPS) detects charged particles scattered at small 
angles and carrying a substantial fraction, xl, of the incoming proton momentum; these parti- 
cles remain in the beam pipe and their trajectory is measured by a system of position sensitive 
silicon micro-strip detectors very close to the proton beam. The detectors are located in six 
stations, SI to S6, placed along the beam hne in the direction of the outgoing protons, at 



^The coordinate system used in this paper has the Z axis pointing in the proton beam direction, hereafter 
referred to as "forward" , the X axis pointing horizontally towards the centre of HERA and the Y axis pointing 
upwards. The polar angle 9 is defined with respect to the Z direction. 
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Z = 23.8 m, 40.3 m, 44.5 m, 63.0 m, 81.2 m and 90.0 m from the interaction point, respec- 
tively. The track deflections induced by the magnets in the proton beam line are used for the 
momentum analysis of the scattered proton. For the present measurement, only the stations 
S4, S5 and S6 were used. With this configuration, for xl close to unity, a resolution of 0.4% on 
the longitudinal momentum and 5 MeV on the transverse momentum has been achieved. The 
transverse momentum resolution is however dominated by the proton beam intrinsic transverse 
momentum spread at the interaction point of ~ 40 MeV in the horizontal plane and ~ 90 MeV 
in the vertical plane. 

The luminosity was determined from the rate of the Bethe-Heitler process, ep e'jp, where 
the photon is measured with a calorimeter (LUMI) located at Z = —107 m in the HERA tunnel 



downstream of the interaction point in the direction of the outgoing positrons [50 



3 Kinematics 

The reactions under study (cf. Fig. |l]) are 

e{k)p{P) e{k')p%V)p{P') and e{k)p{P) e{k')p\V)N{N'), (1) 

where the symbols in parentheses denote the four-momenta of the corresponding particles (or 
particle system, in the case of N). 

The kinematics of the inclusive scattering of unpolarised positrons and protons is described 
by the positron-proton centre-of-mass energy squared, s, and any two of the following variables: 

• Q'^ = —q^ = —(k — k'Y, the negative square of the exchanged photon's four-momentum. 

• y = {q ' ■ P), the fraction of the positron energy transferred to the hadronic final 
state in the rest frame of the initial state proton. 

• W"^ = {q + PY = + '^yi.k ■ P) + — the centre-of-mass energy squared of the 
photon-proton system, where Mp is the proton mass. 

For the exclusive reaction ep — > ep^p tt+tt^) and the proton-dissociative process 

ep ep^N, the following additional variables are used: 

• t = {q — V)"^ = {P — P'Y 1 the four- momentum transfer squared at the photon- p° vertex; 
for the proton-dissociative reaction, t = {q — V)"^ = {P — N')"^. 

• The angle between the p^ production plane (which contains the virtual photon and the 
p°) and the positron scattering plane. 

• The polar and azimuthal angles. Oh and (ph-, of the decay vr^ in the p° helicity frame, 
where the p° is at rest and the polar angle Oh is defined as the angle between the direction 
opposite to that of the outgoing proton and the direction of the vr"*". The azimuthal angle 
(fh is the angle between the decay plane and the p° production plane. 

• xli the fraction of the incoming beam momentum carried by the outgoing proton. 
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For the proton-dissociative reaction, the mass M^v of the diffract ively produced state N 
is relevant. In the present analysis however it was not possible to measure this quantity 
directly and the Mj^ range covered was obtained from Monte Carlo simulations (see 
section 



In the present analysis, events were selected in which the final state positron was scattered 
at an angle too small to be detected in the uranium calorimeter. Thus the angle between the 

production plane and the positron scattering plane was not measured. In such untagged 
photoproduction events, the value ranges from the kinematic minimum Q^j^ = M^y'^/{1 — 
y) ~ 10~^ GeV^, where Me is the electron mass, to the value at which the scattered positron 
is observed in the uranium calorimeter, Ql^ax ~ 4 GeV^, with a median of approximately 
4 X 10~^GeV^. Since the typical is small, the photon-proton centre-of-mass energy can be 
approximated by 

= AEpE.y ~ 2{E, - pzp)E^, (2) 

where Ep, Eg and Ep are the energies of the incoming proton, of the incoming positron and of 
the TT+TT^ system, respectively; the longitudinal momentum of the tt+tt^ system is denoted by 
pzp- Furthermore for = Q^in, t is given by 

t = {q-Vr^ -vlp. (3) 

where p^p is the momentum of the 7r"'"7r~ system transverse to the beam axis. Non-zero values 
of cause t to differ from —"p^p by less than . A multiplicative correction factor determined 
with the Monte Carlo generators discussed in section ^ was applied to the "pip^ distribution to 
account for this effect; the correction was obtained by taking the ratio between the t and p^p 



distributions at the generator level (cf. e.g. |26]). The correction varies from 1.13 at p\ = to 



0.62 at p^ = 0.5 GeV^. The result thus obtained is consistent with that found by using LPS 



tagged events |28|, for which t is measured directly. 



4 Trigger 

ZEUS uses a three- level trigger system [|^, For the present data, the trigger selected events 
from photoproduction of a vector meson decaying into two charged particles without requiring 
that the scattered positron be detected. 

The first-level trigger required an energy deposit of at least 464 MeV in the electromagnetic 
section of RCAL (excluding the towers immediately around the beam pipe) and at least one 
track candidate in the CTD. Events with an energy deposit larger than 1250 MeV in the 
FCAL towers surrounding the beam pipe were rejected in order to suppress proton beam-gas 
events along with a large fraction of photoproduction events. This cut also removes large-MTv 
proton-dissociative events. 

At the second-level trigger, the background was reduced by using the measured time of the 
energy deposits and the summed energies from the calorimeter. 

The full event information was available at the third-level trigger and a simplified recon- 
struction procedure was used. Tighter calorimeter timing cuts as well as algorithms to remove 
cosmic muons were applied. Exactly one reconstructed vertex was demanded, with a Z coor- 
dinate within ±66 cm of the nominal interaction point. Furthermore, the events were required 
to satisfy at least one of the following conditions: 
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1. less than four reconstructed tracks and at least one pair with invariant mass less than 
1.5 GeV (assuming they are pions); 



2. less than six reconstructed tracks with a total invariant mass less than 2.5 GeV (again 
assuming pions). 

Both sets of third-level triggers were prescaled by a factor six. An integrated luminosity of 
2.17 ± 0.03 pb~^ thus yielded approximately 725,000 events. 

No requirements were imposed on the LPS or PRTl at the trigger level. 

For the present analysis, unlike what was done in [EB, the RCAL trigger efficiency at 



the first level was determined |^T], Q using the data rather than a Monte Carlo simulation. A 
sample of two-track events (p" candidates) was used. Since one of the two pions is sufficient 
to trigger the event, the efficiency for RCAL to trigger on a charged pion was evaluated as 
the fraction of events in which the second pion could have satisfied the trigger and in which it 
actually did. This was feasible since, for a subsample of the events, it was possible to uniquely 
determine which of the two pions satisfied the RCAL trigger. The results for the efficiency were 
parametrised as a function of the momentum and polar angle of the pion, separately for positive 
and negative pions. The efficiency was then applied as a multiplicative weight to each event. 
Figure ^ shows the efficiency as a function of the momentum. The uncertainty is dominated 



by statistics. The same parametrisation was also used in ||53 |. 

5 Event selection 

5.1 Event selection for elastic photoproduction 

The following offiine requirements were imposed to select candidates for the reaction ep 



• Exactly two tracks in the CTD from particles of opposite charge, both associated with 
the reconstructed vertex. 

• The coordinates of the reconstructed vertex in the range —0.5 < X < 0.8 cm, —0.8 < 
Y < 0.5 cm and —29 < Z < 38 cm (approximately corresponding to three standard 
deviations of the vertex distribution). 

• Transverse momentum greater than 150 MeV and \ri\ < 2.1 for each of the two tracks, 
thus restricting the data to a region of well understood track reconstruction efficiency. 

• Each CAL cell which is more than 40 cm (in the EMC) or 55 cm (in the HAC) away 
from the extrapolated impact position of either track should not have an energy deposit 
above a given value. The maximum allowed energy deposits varied from 160 to 240 MeV 
depending on the calorimeter part and section. This cut rejects events with additional 
particles, including events with the scattered positron in RCAL. 
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After applying these requirements, the pion mass was assigned to each track and the analysis 
was restricted to events reconstructed in the kinematic region defined by: 

0.55 < M^^ < 1.2 GeV, 

pip < 0.5 GeV^ (4) 
50 < < 100 GeV. 

The restricted range in the two-pion invariant mass Mtttt reduces the contamination from re- 
actions involving other mesons, in particular from production with subsequent — > K'^K~ 
decay and uj tt+tt^tt^ production. The requirement on limits the background from 
proton-dissociative p° production and the selected W range restricts the data to a region of 
well understood acceptance. The final sample contains 79,010 events. 



5.2 Event selection for proton-dissociative yo^ photoproduction 

To select candidates for the reaction ep en^n^N all the criteria discussed for the elastic 
events were applied, except for the cut on the maximum energy deposit in FCAL outside a 
region around the track impact point. In addition, one of the following three requirements was 
imposed: 

• A signal from the PRTl, tagging particles which originate from proton dissociation. A 
signal from the PRTl was defined as a coincidence of signals consistent with that of at 
least a minimum ionising particle from both scintillator counter layers. In addition the 
energy deposit in the FCAL towers around the beam pipe was required to be less than 
1.2 GeV; this was dictated by the trigger condition discussed above. 

This sample contains 2130 events, corresponding to a luminosity of approximately 0.7 pb~^ 
for which the PRTl was operational. 

• An energy deposit in the FCAL towers around the beam pipe between 0.4 GeV and 
1.2 GeV. The lower limit reduces the contribution from calorimeter noise; the upper one 
was again a consequence of the trigger condition. Here also the particles from the proton 
dissociation are tagged. 

A total of 945 events was selected. 

• A proton measured in the LPS carrying a fraction of the incoming beam momentum 
xl < 0.98. As discussed in the xl spectrum measured by the LPS is characterised 
by a narrow peak at x^, ~ 1 from elastic events and a broad distribution for xl < 0.98 
ascribed to proton-dissociative events. The cut xl < 0.98 thus tags the events in which 
the baryon from the proton dissociation is a proton and rejects elastic events, for which 
Xl differs from unity by (Q^ + + |t|)/iy^ [^], i.e. at most 0.2% for photoproduction. 



This sample contains 576 events, corresponding to a luminosity of approximately 0.9 pb ^ 
for which the LPS was operational EH . 



In all cases the region covered is approximately ^10 GeV. This limit is set by the 
requirement at the first-level trigger that less than 1250 MeV be deposited in the FCAL towers 
around the beam pipe (cf. section ^; since M^v could not be measured directly, this limit was 
determined by Monte Carlo studies. 
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6 Monte Carlo generators and acceptance determina- 
tion 



The reaction ep —>■ ep^p was modelled using the EPSOFT generator, developed in the 



framework of HERWIG [Q. The generated Mt^tt, W and t distributions were reweighted so 
as to reproduce the measured distributions after reconstruction. Similarly reweighted were the 
polar and azimuthal angular distributions of the decay pions in the helicity frame. The effective 
W dependence of the 7p cross section was taken as cr oc W^'"^. The t distribution was generated 
as Aexp {—b\t\ + ct^) with 6 = 11 GeV~^ and c = 4 GeV~^. The DIPSI generator was used 



as a cross check of the results obtained with EPSOFT. The LPS acceptance was determined 
using the average of DIPSI and EPSOFT. 

For the simulation of the reaction ep ep'^N, the EPSOFT Monte Carlo was used (in the 
case of the PRTl and FCAL tagged events); for this process, the program is based on the 
assumption that the cross section for the reaction jp —>■ p^N is of the form: 

d (J 1 d(7.yp^pOp I dcpp^pN , d(Tpp^pp\ , > 



dtdM% 2 dt \ dtdMl ' dt 



where the ratio '^"^^^^i^ I '^"^^^^^^ is obtained from fits to pp data 



As a cross check of the results obtained with EPSOFT, the PYTHIA generator p? 



was 



also used (except for the proton-dissociative LPS tagged events for which the acceptance was 
determined with PYTHIA while EPSOFT was used as a cross check). A cross section of the 
form d^ajdtdMl oc e-''l*lF,rf(M^)/M^ is assumed in PYTHIA with 6 = 6o + 2a^' In (W^VM^), 
6o = 2.8 GeV~^ and = 0.25 GeV^^, corresponding to an effective 6^5 GeV~^ in the 
kinematic region for which we present our results. The function ^^^(Mjv) enhances the cross 
section in the low mass resonance region and suppresses the production of very large masses fST 



A fit to the generated Mjv spectrum for 10 < < 200 GeV^ with a function of the type 
gives n = 2.2. The effect of the functions Fsd{M]^) and h = h^M^) on the spectrum 
thus is consistent with the result n = 2.24 ±0.03 obtained for the diffractive dissociation of the 
proton in pp collisions . 



For both EPSOFT and PYTHIA, the value of M^r ranged between Mp+2M^ and a maximum 
fixed by the condition discussed above, Mff/W"^ < 0.1 ||32|, |4^. Although the data extend down 



to Mat = Mp + Mjr, the lack of Monte Carlo events below Mp + 2M^ is not expected to give 
a significant effect The shapes of the two-pion invariant mass distribution and the p° 



decay angular distributions were assumed to be the same as those of the elastic events; this 



assumption is supported by the data, as discussed in section B.l 



The radiation of real photons from the incoming or outgoing positron was not simulated, 
nor were vacuum polarisation loops in the virtual photon; their effects on the cross section were 
estimated to be smaller than 2% 



The generated events were passed through a detailed simulation of the ZEUS detector and 
trigger. They were then subjected to the same reconstruction and analysis programs as the 
data. It was checked that all measured distributions were described well by the simulated 
events. The acceptance in a given bin was then determined as the ratio of the number of 
accepted Monte Carlo events to the number generated in the selected kinematic range. The 
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acceptance, calculated in this manner, accounts for the geometric acceptance, the detector and 
reconstruction efficiencies, the detector resolution and the trigger efficiency. As explained in 
section |^ however, the efficiency of the RCAL trigger was evaluated from the data and then 
applied as a multiplicative weight to each event. 

Figure shows the overall acceptance for elastic events as a function of M^^r, Pxpy cos 9h 
and (fh obtained using EPSOFT. The average acceptance is 15%. The dip in the acceptance 
at the M^^ value corresponding to the p° peak is a consequence of the mass resolution. The 
acceptance as a function of pf.^ and M^, for proton-dissociative events tagged with the PRTl, 
is shown in Fig. ^. As for the elastic case, the acceptance is essentially independent of p'^p. 
While PYTHIA and EPSOFT give consistent results for the shape of the acceptance for the 
PRTl tagged events (and for those tagged in the FCAL or in the LPS with xl < 0.98), the 
normalisation differs by up to a factor of two. 



7 Elastic photoproduction 

7.1 Background to elastic photoproduction 



After the selection cuts described in section 5.1, the main source of background was proton- 



dissociative events in which the mass was small and no particle from the system was 
detected. 



The fraction of proton-dissociative events in the sample selected with the cuts of section E7T 



was determined as follows. Proton-dissociative events were selected with the PRTl (or the 
FCAL, but we shall concentrate on the PRTl tagged sample in the following) as described 



in section p^ . The ratios w of the uncorrected M^^, W, cos9h and iph distributions for the 
proton-dissociative sample (selected with the PRTl) and the sample obtained with the elastic 
cuts (for the period in which the PRTl was operational) were found consistent with being flat, 
as shown in Fig. ^. Since according to both PYTHIA and EPSOFT the requirement of activity 
in the PRTl (or the FCAL) does not affect the shape of the acceptance, this result indicates 
that proton-dissociative and elastic p° photoproduction have the same M^^, W, cos 9h and iph 
distributions. On the contrary, the ratio of the p'^p distributions, also shown in Fig. |^, rises 
with p'^p] since for both reactions the acceptance has the same shape as a function of p'^p (cf. 
Figs. 1^ and H), this indicates a shallower p'^p dependence for the proton-dissociative events. 

The fraction of proton-dissociative events in the total sample was thus taken to depend on p'^p 
only. The p'^p dependence of the background was determined as follows. Let the p'^p distribu- 
tion for the proton-dissociative sample be parametrised as dNdiss/ dp^p = Adiss^W {~K^ssPTp) 
and that for the elastic sample as dNei/dp\p = exp (— fo^^^p^^). Also, let dNpRT-tag/ dp^p 
and dN el-cuts / dp^p indicate the measured p^p distributions for the proton-dissociative sample 
selected with the PRTl and for the sample obtained with the elastic cuts, respectively. Then 
the ratio w{p%p) can be written as: 



w{p. 



_ dNpRT-tag/dpTp 

'^^^ ~ dN,i_cuts/dp^j 



Tp 
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dN diss /dp' 



BeidN^i/dp^Tp + eZsdNd^ss|dplp 

ediss dNdiss/dpTp 

dN^i/dpl^ + {ets^'Veel)dNd^ss/dpl^ 



^ dNdiss/dplp 

dNei/dp'^p + fdissdNdiss/dplp 



Adissexp{-b2ssPTp) 



Aei exp i-b-PPplp) + fdissAdiss exp {-b'^ZPrp 

1 / fdiss 



(7) 



{Aei/fd^ssAd.ss) exp [-{b'^pv - bZs)ph] + 1' 
where Ediss indicates the acceptance for proton-dissociative events to pass the proton-dissociative 



cuts of section |5^ , Sei indicates the acceptance for elastic events to pass the elastic cuts of sec- 
tion p?T| , and Sdi^g"^^^ indicates the acceptance for proton-dissociative events to pass the elastic 
cuts. We introduced the proportionality symbol in eq. (^) to account for the pf^^-independent 
ratio of the acceptance for the proton-dissociative events tagged by the PRTl and that for the 
elastic events satisfying the elastic cuts. The quantity fdiss is the ratio of the acceptance for 
proton-dissociative events passing the elastic cuts and that for elastic events passing the elastic 
cuts; this ratio is taken to be pf^^-independent. 

The difference (6"^^ — &^ff^) was determined by fitting eq. (0) to the data in the range 
< p'^p < 0.5 GeV^ and was found to be 4.8 ± 1.5 (stat.) ± 0.5 (syst.) GeV~^ for the proton- 
dissociative events tagged with the PRTl. The result of the fit is shown as the dashed line 
in Fig. ^. The systematic uncertainty reflects the sensitivity of the result to the limits of 
the fitted range, with the lower limit varied between p'^^ = and 0.075 GeV^ and the upper 
one between 0.3 and 0.5 GeV^. The proton-dissociative events tagged with the FCAL yield 
(b^PP - b'^Z) = 4-1 ± 2.0 (stat.) GeV-2. 

To determine the normalisation of the proton-dissociative background the following proce- 
dure was adopted. As discussed earlier (cf. section ||) the acceptance for the proton-dissociative 
events depends on the Monte Carlo program; hence the proton-dissociative sample was not used. 
Instead we used the sample satisfying the elastic cuts and its subsample ||2^, ^ in which the 
presence of a high momentum {xl > 0.98) proton in the LPS identified elastic events. The re- 
gion 0.075 < p'^p < 0.5 GeV^ was used, where the acceptance of the LPS is well understood ||28[ . 

The function r{p'^p) was introduced, defined as the fraction of proton-dissociative events in 
the elastic sample: 



. 2 N _ dN el-cuts /dp^Tp - dNLPs/dplp 

dNel-cuts/dplp ' 

where dNips/dp'^p is the yield, corrected for the LPS acceptance, for the LPS tagged elastic 
events {xl > 0.98). Using the notation introduced earlier, r{p'^p) can be rewritten as 



"^^^^^ Dexp[-(6'^pp-6rjp|,] + l' 
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where ih""^^ — h^J^^) was taken to be 4.8±1.5 (stat.)±0.5 (syst.) GeV^^, as discussed above, and a 
fit to the data gave D = 7.3lo;9 (stat.) l|j (syst.). In order to correct for the proton-dissociative 
background each event was then weighted by [1 — The resulting integrated fraction of 

proton-dissociative background in the untagged sample is Rdiss = (20 ± 2 (stat.) ^5 (syst.))% 
for pip < 0.5 GeV^. 

In summary, function (^ was fitted to proton-dissociative events tagged with the PRTl to 
determine [b°'PP — fe^ff^). The latter was then used in function (^, which was evaluated using 
the yields for purely elastic events tagged by the LPS {xl > 0.98) and for the events passing 
the elastic cuts. The fit result was used to evaluate the normalisation of the proton-dissociative 
contamination in the sample selected with the elastic cuts and hence the overall contamination 

Rdiss • 

Positron beam-gas and proton beam-gas contaminations were estimated from the unpaired 
bunch event samples to which the selection criteria described earlier were applied. The number 
of events passing the cuts was then scaled by the ratio between the positron (proton) current 
in the paired bunches and the current in the positron (proton) unpaired bunches. The contam- 
ination due to positron-gas interactions was estimated to be (0.6 ± 0.1)%, while that due to 
proton-gas events was found to be ^(0.01 ±0.01)%. The contamination from elastic production 
of uj and (f) mesons (decaying to tt^tt^tt^) was estimated by using simulated events and found 
to be < 0.5% p. 

All subsequent results are shown after subtraction of the contributions from proton-dissociative 
events, beam-gas interactions, u and production. 



7.2 Results for elastic photoproduction 
7.2.1 Cross section determination 

The differential and integrated photoproduction cross sections for the process 7p ir^ir'p 
were obtained from the event yield measured for the reaction ep evr+vr"^. The cross sections 
for these two processes are related by 



ep^eTT^TT p 



a 



dydQ"^ 



'l + il-yf 2{l-y) Q 



2 ' 

min 



y 



27rg2 



I 2(1 ^*p^TT^1T p 



7 p — >7r^7r p 



y 



(10) 



where a is the fine structure constant and a'^^^^ "^{W^Q"^) and ''^'^ ^{W^Q"^) are the 
cross sections for transversely and longitudinally polarised virtual photons, respectively. These 
cross sections are assumed to be 



(11) 
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for transversely polarised photons and 



7 p- 



(12) 



for longitudinally polarised photons, where ^ is a proportionality constant of order unity (cf. 
e.g. m). The results presented in this paper are insensitive to the value of C,"^: varying ^'^ 
between and 1 has negligible effects. 

Substituting the latter two expressions into equation ( ]T0| ) yields: 



' ep— >e7r"'"7r~p 

dydQ'^ 



-p{W{y)), (13) 



which defines the effective photon flux ip{y, Q^). 

From eq. (pTSD, knowing the effective photon flux, it was then possible to determine the 
cross section (T^p^^+,r-p- As an example, the differential cross section dcr^p^^+^-p/dM^jr was 
evaluated in each M^.^ bin of width /S.M„t, as 



'-^^7p— +7r+7r p -^Tr+TT" 



dM^^ A- L-^ ■ AM 



(14) 



TTTT 



where A^tt+tt- is the number of observed events in the bin after background subtraction and 
correction for the RCAL trigger efficiency, L the integrated luminosity and A the overall ac- 
ceptance in the bin excluding the RCAL efficiency. The integral of the effective photon flux 
(/?(?/, Q"^) over the y and ranges covered by the experiment is indicated as In the following, 
for brevity, the subscript 7p vr+vr^p will be dropped. 



7.2.2 Differential cross section da/dM-^T, for the reaction — > tt+tt p 



In Fig. 1^ the differential cross section da/dM^^^^ for the process 7p ti^tt^p is shown in the 
kinematic range |t| < 0.5 GeV^ and 50 < ly < 100 GeV. Figure |^ shows da/dM^^j^ for different t 
bins. The p° resonance shape is skewed, as observed in previous measurements |]l|-[^. This can 
be understood in terms of the interference between resonant production and non-resonant 
Ti~^TT~ production [P3, W^. 



Fits to the points of Fig. ^ were performed in the range 0.55 < 
following parametrisation ||26||: 



< 1.2 GeV with the 



da 
dMI, 



A 



+ B 



+ 1 



P55 



(15) 



where Mp is the nominal p° mass and = TqIp* /pq)^{Mp/Mt^t^) the momentum dependent p° 
width, with Fq the p° width, p* the tt momentum in the tttt rest frame and pg the value of p* 
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Parameter 



value 



stat. error 



To 

A 

B 

Aps 



0.770 GeV 0.002 GeV 

0.146 GeV 0.003 GeV 

-2.75 /ib^/^ 0.04 

1.84 GeV"i/2 0.06 ^b^/^ GeV-^/^ 

0.030 /ib/GeV 0.015 /ib/GeV 



Table 1: Results of the fit to the mass spectrum of Fig. ^ for 50 < W < 100 GeV and 
\t\ < 0.5 GeV^ using expression (|T5|). Only statistical errors are given. 



at the p° nominal mass. The non-resonant amplitude is denoted by B and is assumed to be 
constant and real. The term fps is a first order polynomial of the form fps = Apsil + BpsMj^^^) 
accounting for residual background from the process 7p Xp. The coefficient Bps was 
estimated to be Bps = 1.5 GeV~^ from Monte Carlo studies in which PYTHIA was used 
to simulate the reaction ep eXp. The fitted value of Aps corresponds to an integrated 
contribution typically smaller than 1% of the total, independent of W and p"^. Table gives 



mass 



the parameters of the fit for the spectrum shown in Fig. The fitted values of the p 
and width are consistent with the Particle Data Group tables ||5^ and the /ndf for the fit is 
15.3/21. 

The fits to da/dM-^T^ were repeated in a different way. The mass spectrum was corrected 
for the acceptance excluding the effects of migration. The mass spectrum was then fitted with 
formula (|l^) convoluted with a Gaussian which describes the detector resolution. The width of 
the Gaussian varied between 6 and 14 MeV depending on W. The resulting values for the p° 
mass and width were 771 ± 2 MeV and 159 ± 3 MeV. The difference between this value of the 
width and that given in Table |I| gives an indication of the systematic errors associated to the fit. 
The other parameters of the fit were found to be A = -2.75±0.04 /xb^/^ B = 1.94±0.07 //b^/^ 
GeV-i/2 and Aps = 0.000 ± 0.015 /ib/GeV. The results for the cross section do not change if 
this method is used. A fit including a p^-cu interference term was also performed ||51|; this led 
to a slightly better x^- However, none of the results presented in the following changes if such 
a fit is used. 



The curves shown in Fig. |^ were obtained using a calculation |^ based on Soding's model 
In this case the mass and the width of the p° were fixed to the values given in the Particle Data 
Group tables |E3. The 7r-p total cross section, a free parameter of model [M, was fitted. This 
fit is discussed in detail in section here we only remark that the results of the calculation 
are in good agreement with the data with an average /ndf of 1.0. The non-resonant and 
interference terms are also shown in the figure; the non- resonant contribution varies very little 
with M^-TT, a result consistent with the ansatz, made above, that B is a, constant, as assumed 
in our previous analyses p6|, p8|. 



Fits to the data of Fig. |^ using formula (^) were also carried out, with the mass and the 
width of the fixed to the values given in table |I|. The results for from these fits are 

shown as a function of |t| in the upper plot of Fig. ^: |-B/^| decreases with increasing \t\. 
The quantity |-B/A| is a measure of the ratio of the non-resonant to resonant contributions; 
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its decrease with increasing |t| was already observed in fixed target experiments and can be 



described in the framework of the Soding model p 



Alternatively, the following phenomenological parametrisation proposed by Ross and Stodol- 
sky was used to fit the mass distribution: 



'^'^ --fp- BWp{M^^) ■ {Mp/M^^f + fps, (16) 



dM., 



TTTV 



where EWpIMt^t^) is a relativistic p-wave Breit-Wigner function and the factor [Mp/Mt^t^)^ 
accounts for the skewing of the signal. In this case the fitted values of the p° mass and width 
are 771 ± 2 MeV and 138 ± 3 MeV, respectively; the parameter k is 5.13 ± 0.13. Here again 
the fits were repeated in different |t| bins, keeping the mass and the width of the p° fixed 
to the values just quoted. The parameter k is plotted as a function of \t\ in the lower part of 
Fig. p. The decrease of the amount of skewing with increasing |t| is, in this case, reflected in the 
decrease of k. Our results are in agreement with those found in fixed target photoproduction 
experiments [|l^, [T^ , indicating that skewing of the p° resonant shape depends only weakly, if 



at all, on W . Note that in all t bins the median Q is lower than 10 GeV . The results are 



consistent with the effective expectation of the Soding model p8| , p[] , as the continuous line in 
the lower plot of Fig. || shows (cf. section H). 

Fits using formula ([T5|) were also performed in bins of cos6h and (ph, again with the 
mass and the width of the p° fixed to the values given in table |1|. The ratios |-B/v4| from 
these fits are shown in Fig. \B/A\ appears to be independent of W (as already suggested by 
the comparison with the fixed target data, cf. Fig. ||) as well as of the decay pion polar and 
azimuthal angles in the helicity frame. 



7.2.3 Integrated 7p p^p cross section 



The integrated jp — > p^p cross section for |t| < 0.5 GeV^ was determined in four W bins. 
In each of these bins fits to the mass spectra were performed using equation (|15|); Mp, Tq 
and B/A were fixed to the values given in table |I]. Following refs. 



], the cross section 

was calculated by integrating the resonant contribution obtained from the fit over the range 
2Mt, < Mn-TT < Mp + 5ro. Figure |l^ and Table || show the results. Table ^ gives the results for 
the reaction 7p Ti^n^p over the same mass range; it was obtained by integrating the result 
for the first term in eq. (|l^). 



The systematic uncertainties are dominated by the uncertainties on the acceptance (4-10%), 
the proton-dissociative background (8.5%) and the number of p° signal events, which depends 
on the functional form chosen |5l| to fit the mass spectrum (4%); the parameters Mp, Fq and 
B/A were also varied within their statistical errors (1.5%). The uncertainty on the acceptance 
(4-10%) is W dependent and has two main contributions: the calorimeter trigger efficiency near 
the threshold (4-10%) and the sensitivity of the results to the cuts (4-2%). 

Table H summarises the contributions to the systematic uncertainty. The total systematic 
uncertainty was obtained by summing all contributions in quadrature. 



Figure [TO includes a partial compilation of low energy measurements, as well as the recent 
ZEUS and HI ||2^ results. Also shown are parametrisations based on Regge theory 
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(W) [GeV] 


(T^p^pOp [/ib] 




55 


10.9 ± 0.2 (stat.)Ii 3 


(syst.) 


65 


10.8 ± 0.2 (stat.)ll:? 


(syst.) 


75 


11.4 ± 0.3 (stat.)l{:^ 


(syst.) 


90 


11.7 ± 0.3 (stat.)lji 


(syst.) 



Table 2: Elastic p° photoproduction cross section for \t\ < 0.5 GeV^, integrated over the mass 
range 2M^ < M^-tt < Mp + 5To in four W bins. The results were calculated by integrating the 
resonant contribution obtained from the fit with eq. ([T5|). 



(W) [GeV] 


'^7p— >7r+-7r p [/^b] 


55 


12.2 ± 0.2 (stat.)lli (syst.) 


65 


12.1 ± 0.2 (stat.)l{J (syst.) 


75 


12.8 ± 0.3 (stat.)ll:^ (syst.) 


90 


13.1 ± 0.3 (stat.)j:|| (syst.) 



Table 3: Elastic tt+tt^ photoproduction cross section for \t\ < 0.5 GeV^, integrated over the 
mass range 2Mt, < M^^t, < Mp + 5ro in four W bins. The results were obtained by integrating 
the first term in eq. (0). 



which assume the value of the pomeron intercept found by Donnachie and Landshoff and 



by Cudell et al. |6^, respectively. The W dependence of the data is described satisfactorily by 
both. 

A least squares fit to the present data alone with a function of the type cr^/p^pOp(W) = 
a^p^pOp{Wo){W/Wo)'' gives a^p^pOp{Wo) = 11.2 ± 0.1 (stat.) (syst.) /ib at Wq = 71.7 GeV 
and a = 0.16 ± 0.06 (stat.) Igiis (syst.). The value of a is consistent with the value expected 
for a "soft" pomeron, a ~ 0.22 (see e.g. ||6l|). The systematic uncertainties were determined 
by repeating the fit to the cross section obtained after each systematic check. The differences 
between the values of a^p^pOp{Wo) and a thus found and the nominal value were added in 
quadrature. The dominant contribution to the systematic uncertainty on a is that due to the 
trigger efficiency since its effect in different W bins is not correlated; conversely, the effects of 
the uncertainty of the proton-dissociative background contamination and that of the procedure 
to extract the resonant part of the cross section are the same in all W bins. 



7.2.4 Differential cross section da/d\t 



Figure |TT]a shows the differential cross section d(T/d\t\ for the process 7]? tt+tt^p in the 
kinematic range 0.55 < M^^ < 1.2 GeV and 50 < < 100 GeV. The cross section exhibits 
the exponential fall characteristic of diffractive processes. A fit to the form 



da 

d\t\ 



(17) 
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Contribution from 



Uncertainty 



Luminosity 1.5% 

Acceptance: trigger efficiency 4-10% 

Acceptance: sensitivity to cuts 4-2% 
p-dissociative background subtraction 8.5% 
Background due to elastic u and production 1% 

Procedure to extract the resonant part of the cross section 4.5% 
Radiative corrections 2% 

Total 11-14% 



Table 4: Individual and total contributions to the systematic uncertainty on the integrated 
cross section. 



was performed. The fitted values of 6,r7r and 11.4 ± 0.3 (stat.) 1^:^ (syst.) GeV"^ and 

2.8 ±0.7 (stat.) lj;g (syst.) GeV~^, respectively. The main contribution to the systematic errors 
is the uncertainty of the acceptance. 



In Fig. 12 the slope 67r7r resulting from a fit of equation (^) in different mass bins is shown; in 
this case Ctttt was kept fixed at the value 2.8 GeV""^. The rapid decrease of 6^,^ with increasing 
mass is consistent with the results of earlier measurements (cf. e.g. |T|) and effectively is 
expected in the Soding model ^ as the continuous curve in Fig. |12| shows. The way the 
curve was obtained is discussed in section ^. 



In order to determine da /d\t\ for the resonant process 7p p^p, the mass fits with eq. (p!5| ) 
were carried out in each \t\ bin (with Mp and Fq fixed to the values of table [^) and the 
resonant part of the cross section extracted as a function of |t| and integrated over the range 
2M7r < M^n-TT < Mp + 5Fo. The cross section da/d\t\ thus obtained is plotted in Fig. |Tl|b, where 
the result of the fit with the function 



^ - A p-bp\t\+cpt^ 

d\t\ ~ 



(18) 



is also shown. The parameters of the fit are hp = 10.9 ± 0.3 (stat. 



Cp = 2.7 ± 0.9 (stat. 



+1.9 
-1.7 



^1.0 
-0.5 



(syst.) GeV"2 ^nd 



(syst.) GeV . The larger systematic uncertainty of bp with respect 



to that of ^TrTT reflects the sensitivity to the procedure used to extract the resonant part of the 
cross section. 

Finally the \t\ distribution was studied in three different W bins; the parameter Cp was fixed 
to 2.7 GeV""^. Table ^ and Fig. |13| show the values of the slope bp as a function of W together 
with the other recent results from HERA 



26, 27, BH| and a partial compilation of low energy 



data |, 0, |T9|, (cf. Fig. 9 of ref. [^). 

A fit of the form bp{W) = bp{Wo) + 2a^' \n{W/Wof , with Wq = 71.7 GeV, to the present 
data alone yields a^p' = 0.23 ± 0.15 (stat.) l^o? (syst.) GeV~^. The systematic uncertainty 
was determined by repeating the fit to the b values as modified by the effect of each individual 
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(W) [GeV] 


bp [GeV-2] 


55 


10.6 ± 0.2 (stat.)lo:4 (syst.) 


65 


11.0 ± 0.2 (stat.)lJ:^ (syst.) 


84 


11.1 ± 0.1 (stat.)lJ:[! (syst.) 


Table 5 


bp as a function of W. 



systematic uncertainty 
value were added in quadrature. 



the differences between the values of a^' 



thus found and the nominal 
The present result is consistent with a^' = 0.25 GeV~^ 
obtained from fits to data on soft hadronic processes [^. Such a dependence of bp on W is 
expected to be vahd for > 5-10 GeV pT|. 



7.2.5 Decay angular distributions 

The angular distributions of the decay pions allow one to determine the p° spin density matrix 
elements. They were determined in the helicity frame, where the dependence on 6h and (ph can 
be written as |]63| : 



^ '^'^ ^ [^1 - ^o°o) + ^ fsC' - l) cos^ 9, - v^3fte[r?o^] sin 29, cos ^p. 



ad cos 9hd^h 47r'2V """J 2 

-r^^tiSm^9hCos2ph], (19) 

with Tqq, r°Q and r^ti the p° spin density matrix elements. The element Tqq represents the 
probability that the produced meson has helicity 0; rf^^ is related to the size of the in- 
terference between the helicity non-flip and double flip amplitudes, while 3fJe[r5'Q] is related to 
the interference between the helicity non-flip and single flip amplitudes. If s-channel helicity 
conservation (SCHC) holds, r^i^ and 3fJe[r5'o] should be zero; r^Q should also be small because 
in the kinematic range of the present data the incoming photons are mostly transverse. 



Figure n shows the acceptance corrected 9, and ph distributions for the process jp 



TT^TT p. As discussed above (cf. Fig. their shape is consistent with being the same for 
elastic and proton-dissociative events. 

A two-dimensional least-squares flt of equation (p^ to the acceptance corrected cos9h and 
Ph distributions yields = 0.01 ± 0.01 (stat.) ± 0.02 (syst.), r^t^ = -0.01 ± 0.01 (stat.) ± 
0.01 (syst.) and 3fJe[r°^] = 0.01 ± 0.01 (stat.) ± 0.01 (syst.). The resuh of the flt is shown 
in Fig. |l^. The x^/ndf of the flt is 225/215. A moment analysis gives similar values. The 
systematic uncertainties are dominated by the error of the acceptance. The two-dimensional 
9h, Ph distribution was not corrected for the non- resonant and interference contributions, which 
however appear to have the same cos9h and ph dependence as the resonant process (cf. Fig. ||). 

The present results indicate that in the kinematic range studied the p*^ mesons are produced 
predominantly with helicity ±1. In addition our data are consistent with s-channel helicity 
conservation. 
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The two-dimensional fit described above was repeated in different M^^^ and W bins. The 

in Fig. |T5|; the data do 

M 



results found for Tqq, 3fJe[r°Q] and 
not indicate any dependence on 



,04 



I are plotted as a function of M, 
It should be noted that in some models 



[see e.g. 



for 



finite values of (Q^ > 1 GeV^), rgg is expected to vary with M.^^,. A variation at large values 
of Mtttt, ^tttt ^ 0.9, was observed in fixed target photoproduction experiments [0. The results 
are also independent of W as shown in Fig. 16; for this study the W range was restricted to 
W < 80 GeV since at large W the two-dimensional acceptance in the cos 9h, (fh plane is rapidly 
varying. Here again, the comparison of these results with those obtained for r^^^ and 3fte[r5'Q] by 
the low energy experiments (cf. e.g. |T^) confirms the lack of W dependence for these elements. 
Further investigations, not presented here, show that if SCHC and natural parity exchange in 
the t-channel are assumed, 



then Tqq, 



3?e[rO^] and 

f'l—i appear independent also of t, in the t 



range studied here |51 



8 Proton-dissociative photoproduction 

8.1 Mtttt, W, cosOh, (ph distributions 

As discussed earlier (section [7.1| ), proton-dissociative events selected with the PRTl or the 
FCAL as described in section ^]2| have the same M^^r, W, cos 9^ and iph dependence as the elastic 
events selected by the cuts of section ^]1] (which contain only a contamination Rdiss = 20% from 
proton-dissociative events). This was deduced from the fact that the ratios of the uncorrected 
M^TT, W, cos 9h and (ph distributions for the proton-dissociative sample and the sample obtained 
with the elastic cuts are consistent with being flat, as Fig. | shows. 

This result supports the hypothesis of factorisation of the diffractive vertices. As discussed 
in B^, given the dissociative reaction ha Na and the elastic one ha — ha, with h and a 



hadrons, factorisation implies 



<f(ydiss/dtd{Ml/sha) 
d?aei I dt 



f{s,Mj,/Sha,t), 



(20) 



i.e. at given Sha, M% and t, the ratio of the diffractive dissociation cross section to the elastic 
cross section is a constant independent of hadron a; here Sha indicates the square of the centre- 
of-mass energy of the ha system. 



The results presented in sections ^.2| and |8.3| were obtained for the production of pion pairs 
in the range 0.55 < M^^ < 1.2 GeV and not for the resonant process. This was done because 
of the limited statistics. 



8.2 \t\ distribution 



As discussed in section ITJ] , the and hence the \t\ distribution for proton-dissociative events 
is shallower than for elastic events. The acceptance corrected \t\ distribution for the reaction 
7P — > 7r^7T~N obtained with the PRTl tagged events is shown in Fig. [1^ (solid symbols). 
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The continuous line represents the result of a fit with an exponential function of the form 
j[^-bdiss\t\ jj^ ^]^g range 0.025 < \t\ < 0.5 GeV^ and corresponds to a t-slope bdiss = 5.8 ± 
0.3 (stat.) ± 0.5 (syst.) GeV"^ for the kinematic range 50 < W < 100 GeV and (Mp + M^f < 
M% < O.lW^] the upper limit of {M^^^^ = JOAW^^^ ^ 30 GeV) was chosen following 



refs. [^0|, ^ and corresponds to the region where diffractive interactions dominate. A fit 
with a function of the form A exjp {—b'fHgg\t\ + c^j^s^^) gives b'^^^^ = 6.6 ± 1.1 (stat.) GeV~^ 
and c'^i^, = 1.8 ± 2.4 (stat.) GeV"^ If the analysis is repeated for (Mp + M^f < M% < 
100 GeV^, which is the region covered by the data, the t-slope for a single exponential is 
6.4±0.3 (stat.) ±0.6 (syst.) GeV~^. The dip at low |t| is a consequence of tmin being non-zero at 
large values of Mjy (|tmm| ~ 0.006 GeV'^ for = O.IVV^); it disappears for M^, < 100 GeV^. 



In Fig. O the open circles show the distribution for the LPS tagged events in the kinematic 
range 50 < < 100 GeV and {Mp + M^f < M% < O.lW^. A fit of an exponential function to 
these points yields a slope of 5.8 ± 0.5 (stat.) ± 0.9 (syst.) GeV~^, in agreement with the result 
found with the PRTl tagged events. 

As mentioned earlier, the t distribution determined both for the PRTl and the LPS tagged 
events is for 7p tt^tt^N, not for the resonant process jp — > p^N. From the elastic data 
however one finds that the difference of the t-slopes for the reaction 7p 7i^7i~p and for 
7J9 p^p is ~ 0.5 GeV~^. 

For the result obtained with the PRTl tagged events, the systematic error includes the 
difference of the result obtained with the PRTl and the FCAL tagged events, as well as the 
sensitivity to the Monte Carlo model used (EPSOFT vs PYTHIA) and to the shape of the 
generated da/dM% oc (I/Mat)" spectrum (n was varied in the range 2.0 < n < 2.4). For the 
result obtained with the LPS tagged events, the sensitivity to the selection cuts and the fitted 
\t\ range was also included. 

The t-slope in proton-dissociative p^ photoproduction is thus smaller than that for the elastic 
process by about 5 GeV~^. This is consistent with the results found for virtual photons 



and with those obtained for hadron-hadron collisions |^2[-||33]; it is also in agreement with 
theoretical estimates (cf. e.g. 



8.3 Ratio of the elastic to the proton-dissociative photoproduc- 
tion cross sections 

Because of the large discrepancy in the normalisation of the acceptance obtained with PYTHIA 
and EPSOFT for proton-dissociative events, the cross section for the process 7p — > p^N was 
not determined directly from the PRTl, FCAL or low xl LPS tagged events. We instead 
determined the ratio of the elastic to the proton-dissociative p° photoproduction cross sections, 
a^p^pOp/a^p^pO^, using the ratio Rdiss found by means of the LPS tagged events with xl > 0.98 
(cf. section [7. 1| ) . 

The elastic yield was calculated as N{1 — Rdiss) / ^eu where N is the number of events passing 
the selection criteria presented in section |5.1| and Rdiss is the fraction of proton-dissociative 
events in this sample (see section [7.1| ). The acceptance Sei is the one determined with EPSOFT 
for elastic events. The proton-dissociative yield was determined from N Rdiss / ^'dilT*'^ ) where 
^diss^^'^ is the acceptance for proton-dissociative events when the criteria used to select the 
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elastic events are applied (section Note that for e^J-^""^' EPSOFT and PYTHIA agree in 
shape and normalisation. For the kinematic range 50 < < 100 GeV, \t\ < 0.5 GeV^ and 
(Mp + M^)2 <Ml< Q.lW^ one obtains 



■el/ p— diss 



Rdi 



diss ^diss 



_el—cuts 



cr. 



Rdi. 



'diss 



£el 



2.0 ±0.2 (stat.) ±0.7 (syst.). 



(21) 



The quoted error is given by the uncertainties on Rdiss and on the acceptance. The result 
was obtained assuming a mass dependence of the type cPa/dM^^ oc l/M^ with n = 2.24 as 
measured by CDF Varying n by ±0.2 induces a change of ±0.3 in (J-yp^pOp/ajp 
is not included in the quoted systematic uncertainty. 



pOjv; this 



The present result is consistent with that found for pp collisions at ISR p3, R, 



2.08 ± 0.13 at a centre-of-mass energy of 53 GeV for < 0.05spp 



'^el/p—diss 



is the proton-proton 

centre-of-mass energy) and 1.69 ± 0.11 for < 0.1 Spp. It is also consistent with the result 
found by the HI Collaboration (T^p^pOp/cr^p^^pOjv = 1-54 ± 0.26 (stat.) ±0.31 (syst.) for 
7 < g2 < 36 GeV2, 60<W < 180 GeV and Mjj < O.ObW^. The ZEUS and HI data together 
indicate that the ratio (Jyp^pOp/cr^p^pQN is not a strong function of the photon virtuality. Our 
result in conjunction with the pp data and the ep results at non-zero supports the hypothesis 
of factorisation. 



9 A model dependent derivation of the pion-proton 
cross section 

As discussed earlier, the measured cross section for the process '-/p Tr+vr^p includes the 
contributions of resonant vr+vr^ production, non-resonant vr+vr^ production and their 

interference. Non-resonant tt^tt^ production can be described by the photon fluctuating into a 
virtual tt+tt^ pair with one or both pions scattering elastically off the proton. The amplitude 
for this process can thus be written in terms of the rep total cross section a^rp- We have extracted 
this cross section in the framework of a recent calculation [^, based on Soding's approach ||38 
in which aj^p is one of the parameters. 



The total 7cp cross section was determined by fltting the calculation of ref. |^T| to the Mt^^ 
distribution of Fig. |^. The flt gives a^p = 31 ±2 (stat.) ±3 (syst.) mb at an average pion-proton 



centre-of-mass energy y^s^ ~ y {W'^)/2 ~ 50 GeV. The systematic error reflects the systematic 
uncertainty of the data. An additional uncertainty of approximately 15% was evaluated by 
repeating the flt with different values of the parameters of the model. The value of x^/ndf is 
23.4/23. 

Our result is consistent with the extrapolation of the flts by Donnachie and Landshoff 
which give a-^p = 26.6 mb at w^s^ = 50 GeV. 



The predictions of the calculation [O] using the fltted value of shown in Figs. 0, |] 

and[l^. Both the decrease of the skewing with increasing \t\ and the variation of the t-slope with 
M^n-TT are well described. To obtain the curves shown in Figs. B| and 112], events were generated 



with a Monte Carlo program based on and were binned as a function of M^tt and t. The 
flts performed to the Mtttt spectra in the data for different t bins (using eq. (|T^)) and to the t 
distributions for different M^^t^ bins were repeated for the generated events. 
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10 Summary and conclusions 



We have presented a high statistics study of p° photoproduction for 50 < ly < 100 GeV and 
\t\ < 0.5 GeV^. With respect to previous analyses at HERA, the present one features larger 
statistics and reduced systematic uncertainties. The main novel results can be summarised as 
follows: 



The TT+TT invariant mass spectrum is skewed and the amount of skewing decreases with 
increasing consistent with the results from fixed target experiments. 

The cross section for resonant production, 'jp p^p, is 11.2 ±0.1 (stat.) tli (syst.) /ib 
at (W) = 71.7 GeV. It increases slowly with W, exhibiting a power-like behaviour of the 
type W°- with a = 0.16 ± 0.06 (stat.) l:o:i5 (syst.), consistent with a ~ 0.22, the value 
expected for a "soft" pomeron. 

The t distribution for the reaction 7p Ti^n^p is well described by an exponential of 
the form A^^^^exp (— + CT.T.t'^). The slope Bt,^ decreases rapidly with increasing Mj^^, 
again consistent with the results from fixed target experiments. 

The t dependence of the cross section of the reaction 'jp —>■ p^p can also be described by a 
function of the type Apexp {-bp\t\ + Cpt^), with bp = 10.9 ± 0.3 (stat.) (syst.) GeV^ 
and Cp = 2.7 ± 0.9 (stat.) tj;? (syst.) GeV'^ 

A fit with the function bp{W) = bp{Wo)+2a^' \n {W/Woy yields a^' = 0.23 ±0.15 (stat.) 

+0.10/ 
-0.071 



'^■^'^.(syst.) GeV ^, consistent with results from elastic hadron-hadron scattering. 



The spin density matrix elements Tqq, 3?e[r°o] and r^'^^ were obtained from the angular 
distributions of the decay pions in the helicity frame; their values are consistent with 
s-channel helicity conservation. No dependence on M-^t^ or W is observed. 

Proton-dissociative p" photoproduction, in which the proton diffractively dissociates into 
a system with mass M^v < 10 GeV, exhibits dependences on M^^^, W, cos9h and iph 
consistent within errors with those of the elastic process. The slope of the t distribution 
is smaller than in the elastic reaction and for 0.55 < M^-vr < 1.2 GeV and {Mp ± M^r)^ < 
Mfj < O.IW^ is measured to be 6 = 5.8 ± 0.3 (stat.) ± 0.5 (syst.) GeV"^ using the PRTl 
tagged events. In this kinematic region, the ratio of the elastic to proton-dissociative 
cross sections is 2.0 ± 0.2 (stat.) ± 0.7 (syst.). 

A model calculation ^l| based on the Soding approach ||38| was fitted to the M^^ spectrum 



for the reaction 7p n^-K^p. The fit yielded a^^p = 31 ± 2 (stat.) ± 3 (syst.) mb at 
an average pion-proton centre-of-mass energy of approximately 50 GeV, consistent with 
the predictions of fits to fixed target np data based on the "soft" pomeron. The model 
dependent uncertainty was estimated to be approximately 15%. 



In p" photoproduction the photon thus appears to behave like a vector meson. The W and t 
dependences of the cross section are those expected for elastic hadron-hadron scattering and the 
object mediating the interaction appears to be the same pomeron that dominates the hadron- 
hadron total cross section. The comparison of the elastic and the proton-dissociative reactions 
suggests that the coupling of the pomeron to the photon is independent of that to the proton. 
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as expected on the basis of factorisation. Our results indicate that the p° not only carries the 

quantum numbers of the photon, but also its helicity in the s-channel system is equal to that 
of the photon. The skew of the mass shape and its t dependence can also be understood in 
terms of soft hadron-hadron interactions and simple quantum-mechanical interference between 
resonant and non-resonant production of pion pairs. 
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Figure 1: Elastic or proton-dissociative p production in ep collisions. 
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Figure 2: RCAL trigger efSciency as a function of the pion momentum P^^. The full symbols 
refer to positive pions and the open ones to negative pions. Only statistical errors are shown. 
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Figure 3: Overall acceptance e for elastic p photoproduction, ep — >■ ep p, as a function of 
(a) Mtttt, (b) (c) (d) cosOh and (e) obtained with the EPSOFT generator. Only 
statistical errors are shown. 
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Figure 4: The overall acceptance as a function of (a) pf.^ and (b) M^gen for proton-dissociative 
events, ep — > ep^N, tagged with the PRTl (see section |5^ ); Mffgen indicates the generated 
value of Mat. Only statistical errors are shown. 



29 



ZEUS 1994 



w 



0.1 



0.05- 



w 



0.1 



0.05- 



0.5 0.75 1 

M„, (GeV) 

-t^+i\ 



-1 

W 0.2 — 
0.15- 
0.1 - 
0.05 



1 
COS 0. 




80 100 
W (GeV) 



w 
0.1 



0.05 







2 4 6 



• m ^ 



0.1 



0.2 



0.3 



0.4 0.5 

Ptp (GeV-) 



Figure 5: The ratios of the uncorrected M^-tt, cos6'ft, y^/i and p^^ distributions for the 
proton-dissociative sample (PRT tagged) and the sample obtained with the elastic cuts. Only 
statistical errors are shown. The dashed line is the result of the fit with equation (|^). 
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Figure 6: The differential cross section da/dM^^T^ for the elastic reaction 7p — >■ tt^tt'p in the 
kinematic region 50 < ly < 100 GeV and \t\ < 0.5 GeV^. The points represent the ZEUS data 
and the curves indicate the result of the fit to the data using expression (plsl). The dashed curve 
represents the resonant contribution, the dot-dashed curve the non-resonant contribution and 
the dotted curve the contribution of the interference term. The continuous curve is the sum. 
Only statistical errors are shown. 
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Figure 7: The differential cross sections da/dM^^T^ for the elastic reaction 7p — tt^tt^p in the 
range 50 < < 100 GeV for different \t\ bins. The points represent the ZEUS data and the 
curves indicate the results of the fit to the data based on the model of ref. |^ (cf. section 
The dashed curve represents the resonant contribution, the dot-dashed curve the non-resonant 
contribution and the dotted curve the contribution of the interference term. The continuous 
curve is the sum. Only statistical errors are shown. 
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Figure 8: The ratio |-B/A| (upper plot) and the parameter k (lower plot) as a function of |t| 
obtained by fitting eq. ( 1T5D and eq. (p!6D, respectively, to the points of Fig. ^ for the elastic 
reaction jp tt^tt^p in the region 50 < W < 100 GeV. The solid points represent the ZEUS 
measurements. The inner error bars indicate the statistical uncertainty, the outer ones the 
statistical and systematic uncertainties added in quadrature. The results of the fixed target 

13 (labelled "Ballam et al." 



experiments 
also shown. 



T|] and 



and "Gladding et al.", respectively) are 



The continuous line indicates the effective expectation of the Soding model as 
implemented in the calculation by Ryskin and Shabelski (cf. section ^). 
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Figure 9: The ratio as a function of W, cos 6h and for the elastic reaction tt^tt p 

in the kinematic range 0.55 < M^^ < 1.2 GeV, 50 < W < 100 GeV and \t\ < 0.5 GeV^. 
Statistical errors only are shown. The continuous lines indicate the average value of 
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Figure 10: The integrated cross section a. 



IP- 



function of the centre-of-mass energy W. 



The ZEUS results are given for the range 2Mn- < M^^r < Mp + 5Tq 



results from HERA |2g, ^ and a compilation of low energy data [Il4["119l) H 



t| < 0.5 GeV^. The other 
are also 



shown. The continuous and dashed line are parametrisations |B2[ based on Regge theory which 
assume the value of the pomeron intercept found by Donnachie and Landshoff and by 
6^ , respectively. The band corresponds to the uncertainty in the determination 

The error bars of the ZEUS points indicate the sum of 
statistical and systematic uncertainties in quadrature. For the points at the same value of W, 
the error bars have been offset. 



Cudell et al. 
of the pomeron intercept of ref. |E2 



35 



ZEUS 1994 




Figure 11: (a) The differential cross section da/d\t\ for the process 7p — > tt^tt'p in the kine- 
matic range 0.55 < M^^ < 1.2 GeV and 50 < < 100 GeV. (b) The differential cross 
section da /d\t\ for the process 7p p^p in the kinematic range 2M^ < M^-^ < Mp + 5ro and 
50 < W < 100 GeV. The continuous lines in (a) and (b) represent the results of the fits with 
the functional forms ([T7|) and (0), respectively. The error bars represent only the statistical 
uncertainties and are smaller than the size of the symbols. 
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Figure 12: The slope ^TrTr resulting from a fit of equation (0) to the t distribution for the 
reaction '-/p — > ir'^n'p in different mass bins. The kinematic range is 50 < W < 100 GeV 
and \t\ < 0.5 GeV^. The continuous line indicates the effective expectation of the Soding 
model as implemented in the calculation by Ryskin and Shabelski (cf. section y). The 
inner bars indicate the statistical uncertainty and the outer ones the statistical and systematic 
uncertainties summed in quadrature. 
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Figure 13: The slope bp for the elastic reaction p p in the kinematic region 

50 < W < 100 GeV and \t\ < 0.5 GeV^ as a function of W together with the other re- 
cent results from HERA ^ and a compilation of low energy data p5|, |1^, |19|, |17 



The continuous line shows the result of the fit discussed in the text; the extrapolation of the fit 
to the low W region is indicated by the dashed line. The error bars of the HERA data indicate 
the statistical and systematic uncertainties summed in quadrature. For the points at the same 
value of W, the error bars have been offset. 
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Figure 14: The differential distributions {1 / a) {da /d cos Oh) and {1 /a) {da/ dip h) for tlie reaction 
7P — TT"*"??"!) in tlie kinematic region 50 < < 100 GeV and |t| < 0.5 GeV^. Tlie continuous 
lines represent the results of the fit discussed in the text. Only statistical errors are shown. 
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Figure 15: The results for Tqq, 3?e[r°Q] and rji-,^ as a function of M^ttt for the reaction 7p — tt^tt p 
in the kinematic range 0.55 < M^^ < 1.2 GeV, 50 < W < 100 GeV and \t\ < 0.5 GeV^. The 
inner bars indicate the statistical uncertainty and the outer ones the statistical and systematic 
uncertainties summed in quadrature. 
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Figure 16: The results for Tqq, Ke[r°Q] and rf^^ as a function of W for the reaction 7p — > tt+tt p 
in the kinematic range 0.55 < M^^ < 1.2 GeV, 50 < < 80 GeV and \t\ < 0.5 GeV^. The 
inner bars indicate the statistical uncertainty and the outer ones the statistical and systematic 
uncertainties summed in quadrature. 
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Figure 17: t distribution for the reaction 7p tt+tt N tagged with the PRTl (full sym- 
bols) and with the LPS {x^ < 0.98, open symbols) in the region 0.55 < M^^r < 1-2 GeV, 
50 <W < 100 GeV, |t| < 0.5 GeV^ and (Mp + M^f < < OIW"^. The dip at low \t\ is 
a consequence of tmij; being non-zero at large values of Afjv- Only statistical errors are shown. 
The line shows the result of the fit to the PRTl points described in the text. 
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